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Abstract 



The running of the effective electromagnetic couphng is for many electroweak 
observables the dominant correction. It plays an important role for deriving 
constraints on the Standard Model in the context of electroweak precision 
measurements. We compute the four-loop QED corrections to the running of 
the effective electromagnetic coupling and perform a numerical evaluation of 
the different gauge invariant subsets. 
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1 Introduction 



The input parameters of Quantum electrodynamics (QED), like the fine-structure con- 
stant or the masses of the charged leptons are known with high precision. For example, 
the fine-structure constant has been measured at the level of ~ 0.32 ppb [1] or the mass 
of the electron in units of MeV is known at the level of ~ 22 ppb [1]. Also from theory 
side physical quantities in QED can be computed to high accuracy in perturbation theory, 
since the electromagnetic coupling constant is a small parameter. 

QED corrections play an important role for electroweak precision observables where 
the effective electromagnetic coupling at the Z-boson mass scale enters. The change 
from the fine-structure constant a, defined at zero momentum transfer, to the effective 
electromagnetic coupling at the Z-boson mass scale is mediated by the photon vacuum 
polarization function. At the Z-boson mass scale the effective electromagnetic coupling 
is governed by the parameter Aa, which plays an important role in constraining the 
Standard Model (SM) in the context of electroweak precision measurements, so that its 
precise theoretical understanding is an important task. 

The parameter Aa receives leptonic and hadronic contributions. The latter escape a 
perturbative treatment and are determined using cross-section measurements of electron 
positron annihilation into hadrons, see e.g. Refs. [2Hn]- The hadronic contributions Aohad 
constitute currently the dominant uncertainty to Aa. The leptonic contributions Aaiep at 
one- and two-loop order in QED are known from the results for the vacuum polarization 
function of Refs. [I0l[ll]. The three-loop order has been computed in Ref. |12j in terms 
of expansions in the small mass ratio M|/M|, where Mi stands for the masses of the SM 
charged leptons and Mz is the Z-boson mass. The dominant contribution to Aaiep arises 
form large logarithms in this small mass ratio, so that it is often sufficient to perform the 
expansion up to the constant term, including it, and to neglect power suppressed terms of 
higher order in the mass ratio M|/M|. Within this work we will extend the calculation 
of Aaiep to four- loop order in QED, also in form of small lepton mass expansions up to 
the constant term. For the lower orders in perturbation theory we include the power 
suppressed terms. Next to the effective electromagnetic coupling we will also present an 
numerical evaluation of Aaiep, which allows to study the size of the different contributions. 
As a by-product of our calculation we will also provide the result for the relation between 
the MS and on-shell mass from the contribution due to heavy fermion loops at three-loop 
order in QED. 

In the next Section [2] we start with some generalities and the definition of our notation. 
In Section |3] we give an outline of our calculation. The new results for the vacuum polariza- 
tion function in the on-shell scheme at four- loop order are presented in Section H] including 
a numerical evaluation of Actiep. Finally we close with our summary and conclusions in 
Section [5l In the Appendix we provide known results as supplementary information. 
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2 Generalities and notation 



We consider the photon vacuum polarization in QED with all one particle irreducible 
insertions into the photon propagator. The vacuum polarization tensor is given by 

W^{q,M) = z I dxe''^^0\Tf{x)fmO), (1) 

where j'^i^x) is the electromagnetic current. It can be decomposed with respect to its 
Lorentz structure into two terms 

w{q,M) = (-g V + gV) n (g^ M) + gvnL (g^ M) , (2) 

where 11 (g^, M) and 11^ (g^, M) are the transversal and longitudinal parts of the vacuum 
polarization. The longitudinal part vanishes due to the Ward-Takahashi identities. In the 
following the symbols Mg, M^, denote the on-shell masses of the electron(e), muon(/i) 
and tauon(r) with the mass hierarchy <K <<c M^-. The polarization function 
n(g^,M) depends on the lepton masses M = {Me,Mn,Mr} as well as on the external 
Minkowskian momentum q, where we consider here > 0. It is related to the effective 
electromagnetic coupling by 

which describes the running of the electromagnetic coupling to higher energy scales. The 
running from g^ = to the scale of the electroweak Z-boson mass, g^ = M|, is mediated by 
the parameter Aa, where we consider within this work the purely leptonic contributions, 

Aaiep = Aaiep(M|) = -ReHoslg' = M|) , (4) 

from the Standard Model leptons which we defined above. The symbol Re denotes the 
real part of a complex quantity and the label os specifies that the vacuum polarization 
function has been renormalized in the on-shell scheme. In the on-shell scheme holds 

nos(g' = 0) = 0. 

In order to determine Aaiep we need to compute the vacuum polarization function in the 
on-shell scheme. We define its perturbative expansion as 

^os(g^M) = £5^(^)'^g)(g^M), (5) 

fc=0 

where each term on the r.h.s. of Eq. ([5]) stands for the one-loop(A; = 0), two-loop(fc = 1), 
three-loop(fc = 2) and four-loop(fc = 3) contribution. The vacuum polarization function in 
the on-shell scheme at one- and two- loop order is known since long p!0|[TT] . The one-loop 
result reads 
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with (5f = I - Xi, Xi = {2Mif/q^ and L^m, = log(-gV^?)- The elhpsis in Eq. ^ 
represent the contributions from higher orders in the expansion in Mf / for an external 
momentum which is much larger than the threshold for lepton pair production. The 
index i stands for the flavor of the lepton which is connected to the external photons. 
This notation will be kept also in the subsequent equations. The analytic continuation of 
the logarithm is determined by the ie-prescription. 
The two-loop result is given by 
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The diagrams which contribute to the one- and two-loop vacuum polarization function 
are shown in Fig. [H 




(a) (bi) (b^) (ba) 

Figure 1: Diagrams contributing to the one-loop and two- loop vacuum polarization func- 
tion. The solid lines denote charged leptons (e, /i, r); the wavy lines denote photons. 

The corresponding result in QCD, where quarks instead of leptons are running in the 
fermion loop, can easily be deduce up to two-loop order from Eqs. ([6]) and ([8]) by multi- 
plying them with the appropriate color factors. 

The three-loop contribution has been determined in terms of expansions in the low 
(g^ — )■ 0) and the high energy limit in Refs. |13fET| . The result in the on-shell scheme has 
been computed in Ref. ^2], where it has been decomposed into four contributions 

^S(g^M) = 5^ (^J)(g^M,) + ^g)(g^M,)) 

+ nfH?',M„M,) + Y (g^M„M,). (9) 

i=^^T i=e,fi 
3+^ 
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The first term H)^ (g^,Mj) is the quenched contribution which originates from diagrams 
without any insertion of an internal, closed fermion loop. In contrast, the terms Il^p\q'^, Mi), 
np^(g^. Mi, Mj) and 11^^^ (g^. Mi, Mj) arise from diagrams with internal lepton loop inser- 
tions. The contribution with the subscript F arises from diagrams with an internal lepton 
loop of the same flavor as the outer loop which is connected to the external photons; the 
contribution with the subscript h comes from diagrams with a heavier lepton inserted in 
the inner loop than in the outer one, and the subscript / denotes the contribution where 
the mass of the lepton in the inner loop is lighter than the external one. An example- 
diagram for each of these four different kind of contributions is shown in Fig. O 




Figure 2: Examples for the different diagram classes at three-loop order. For each class 
one typical representative is shown. Loops of particles which are represented by a dashed 
line denote particles which have a lower mass than the ones which are represented by a 
solid line. 

The four different contributions on the r.h.s. of Eq. ([9]) are renormalized in the on-shell 
scheme and we omit a label indicating it for brevity. The index i labels again the ex- 
ternal leptons which connect to the external photons in a diagram, whereas the index j 
enumerates the internal leptons. The three-loop on-shell results for the different terms 
on the r.h.s. of Eq. were determined in terms of expansions in the limit » {2Mf 
in Ref. [12] and are shown for completeness in Appendix lA.H where we have in addition 
also added mass corrections of higher orders in the small ratio M^/g^. For n^'*(g^,Mj) 
and ng^(g2,Mi) they can be obtained from the MS results in Refs. [T6l|l71[l9ll2l] after 
converting them from the MS to the on-shell scheme. 

The four- loop result can be decomposed, similarly to Eq. (j9]), into several gauge in- 
variant contributions according to the properties and number of the inserted lepton loops 

ngkg',M) = Yl [^S^(9^M,) + ^g].(g^M,) + ^g)(g^M,)" 

i=e,fi,T 

+ Yl [n?(^'' M,) + 4t (g^ M„ M,) + M,)_ 

i=e,ii 
i=ti,T 
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The property which characterizes each individual term on the r.h.s. of Eq. ffTOl) is defined 
in Fig.[3]by the shown example-diagrams. Singlet diagrams, in which the external photons 
are not connected to the same fermion loop did not yet contribute at three-loop order 
due to Furry's theorem [22] and survive for the first time at four-loop order. They are 
given by the terms n^^](g^, M,) and nf^^ [q^-, max{Mfc, M/}), where the subscript d stands 
for 'diagonal' which denotes diagrams with the same lepton fiavor in the two fermion 
loops, whereas the subscript nd stands for 'non-diagonal' and denotes diagrams with 
different lepton fiavors in the two fermion loops. The diagrams with three different lepton 
flavors which correspond to the three terms li^p with j = {e/i, er, /ir} of Eq. f lTOj) . also 
arise for the flrst time at four-loop order. As in the previous Eq. the subscripts 
{FF,F,h,Fh,hh,l,Fl,ll} of the non-singlet contributions in Eq. f lTU]) refer always to the 
internal leptons. Apart from that we show only the function arguments of the lepton 
masses on which these expanded results will depend to the considered expansion depth. In 
general, without expansion, the functions may depend on more mass arguments than those 
which are given here. The on-shell results for Uf{q'^,Mi), IlPp{q\Mi), Ilf{q'^,Mi) and 
n^'^](g^. Mi) have already been determined in Refs. [23l|2l] in the high energy expansion. 
They are given for completeness in Appendix lA. 21 The calculation of the terms in Eq. ( ITOl) 
which are still unknown in order to describe the complete lepton mass hierarchy will be 
discussed in the next section. 



In order to compute the leading four-loop QED corrections to Aaiep we will determine 
the on-shell vacuum polarization function contribution nQg(g^ = M|,M) of Eq. ( fTOl) 
in the small lepton mass expansion up to power corrections of the order M^/M|. For 
this purpose we need to know the MS vacuum polarization function on the one hand for 
massless leptons and on the other hand we need to know it in the low energy limit at 
q^ = for the renormalization procedure in the on-shell scheme as we will see in the 
following. 

The vacuum polarization function in the massless limit, 11 [q\m = 0), is known 
analytically in the MS scheme up to four-loop order (A; = 3) [T3tfTl[2nl[23ti27]. In this 
limit it does not depend on the lepton masses but only on the external momentum. The 
perturbative expansion of the vacuum polarization function in the MS scheme is deflned 
here in complete analogy to the on-shell case shown in Eq. (|5]), only that all quantities 
are renormalized in the MS scheme. 

The non-singlet contributions of the vacuum polarization function in the MS scheme 
at four-loop order in QCD and at = were computed in Ref. [28]. The QED result 




3 Calculation 
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was presented in Refs. [23|[2^. Higher moments of the low energy expansion in /rrf' 
were determined in Refs. [281 - IM] at four-loop order. In these works the contributions 
from quenched diagrams were determined as well as the contributions from diagrams 
with inserted fermion loops, where the internal fermions were massless or had the same 
mass as the external one. In order to perform the determination of Aaiep in the previously 
described expansion at four-loop order in QED one needs in addition still to compute the 
contributions at = from diagrams where the internal fermion loops have a larger 
mass than the external one as well as the non-diagonal singlet contributions. In Fig. [3] we 
show one four-loop example- diagram for each arising gauge invariant subset of diagrams. 
These gauge invariant subsets are given by the different terms in Eq. ( ITOl) . 
The two ingredients, the vacuum polarization function in the MS scheme in the massless 
limit and at = 0, are then combined to construct the vacuum polarization function in 
the on-shell scheme 





+ 
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k=0 

(11) 

where the 'bar' on the r.h.s. indicates that a quantity has been renormalized in the MS 
scheme, e.g. a is the fine-structure constant in the MS scheme. Masses in the MS scheme 
are denoted by small letters fn, whereas masses in the on-shell scheme are given by capital 
letters M. The conversion function CfnM{Caa) converts the lepton mass(fine-structure 
constant) from the on-shell to the MS scheme or vice versa. 

In order to compute 11 (g^ = 0, m) at four-loop order in QED we start first to generate 
the required diagrams with the program QGRAF [35j. The expansion around = is 
performed with the programs q2e and exp [36VI38] which also map the resulting loop 
integrals onto a proper notation, which allows us to perform further simplifications of the 
resulting loop integrals in a straightforward way. In the next step the integrals are reduced 
to a small set of master integrals with the traditional integration-by-parts(IBP) method 
in combination with Laporta's algorithm [39lll0]. For the algebraic manipulations we use 
FORM [¥Tlfl3] and the rational functions in the space-time dimension d, which arise while 
solving the linear system of IBP equations, are simplified with FERMAT [H] . At the end of 
this procedure all four-loop integrals can be expressed in terms of 12 master integrals from 
which three are factorized. They are known since long in the literature [IHlSSj; also other 
authors have contributed to individual master integrals or specific orders in the e = 2—d/2 

expansion with analytic results [29P71 - I53] which allow to obtain 11 (g^ = 0, m) in Eq. f lTTj) 
completely analytically. 

To check our calculation we have kept the dependence on the gauge parameter of 
the photon propagator and have verified its cancellation in the final result. We have 
also checked the Ward-Takahashi identity and verified that the longitudinal part of the 
vacuum polarization function vanishes. 

For the conversion of the vacuum polarization function to the on-shell scheme in 
Eq. ( fTTl) the conversion functions CmM for the lepton mass and Caa for the fine-structure 
constant are needed. 
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Figure 3: For each of the different terms in Eq. f lTOj) is shown here one typical four- loop 
diagram as representative of the whole class. The wavy lines denote photons; the lines 
with an arrow are leptons. Dashed lines represent particles which have a smaller mass 
than particles which are given by solid lines. Dotted lines denote particles which have 
a smaller mass than particles which are represented by dashed lines. The singlet type 
diagrams (11^ ^, 11^ contribute for the first time at four-loop order in QED. 

Let us start with the conversion of the lepton masses. For the case of the mass of the 
r-lepton the QED MS-on-shell relation is known to three-loop order; it can be deduced 
from the QCD results of Refs. [5H460] . We decompose the different orders in perturbation 
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theory as follows 



m£ = Mi 



(5)%f' + o(.*) 



(12) 



The coefficient functions cf^ (/c = 1, 2, 3) depend in general on the renormalization scale /i 
and on the lepton masses, where we have omitted these function arguments for simplicity. 
The one- loop coefficient is known since long [54j and given in Eq. ( I47l) . Starting from 
two- loop order also internal lepton loops can arise, like shown in Fig. HI In general, the 

(2) 

coefficient c\ can be decomposed into gauge invariant contributions which arise from 
quenched {A) diagrams, contributions from diagrams with light (/) or heavy (h) inserted 
lepton loops and diagrams where the internal lepton loop has the same mass as the 
external one {F). 

The tauon is the heaviest lepton of the SM, so that only light internal leptons or the 
tauon itself can contribute in its conversion relation of Eq. fll2p . At two-loop order one 
obtains 

(M.) = c^^\m,) + 2cf\M,) + cf{M,) . (13) 

The individual terms on the r.h.s. of Eq. (fT3|) are known [55] and are given in Appendix [Cl 
At three-loop order the same decomposition with respect to the number of inserted 
closed lepton loops can be made 

c(3)(M.) = cf{Mr) + 2cf )(M.) + cfiMr) + 2cg?(M.) + 4cf (M.) + cg],(M.),(14) 

where the last three terms with two letters in the subscript arise from diagrams with two 
lepton loop insertions. The coefficients on the r.h.s. of Eq. ( |T4l) are again known [56l 
[57l[59] - [6T] and given for completeness in Appendix O In particular the complete effects 
from a second virtual massive fermion in the inner loops were determined at two- and 
three-loop order in Refs. [55l[60]. Since our final result for Aaiep at four- loop order will 
be in an expansion in M^/M| up to the constant order, we do not consider higher order 
virtual mass effects in Eqs. (IT^ and (fn|) as well as in the following Eqs. (JISD, (ITBi) . (HM 
and (|20|). 

For the conversion relation of the muon internal heavy r-lepton loops can arise starting 
from two-loop order. These internal r-lepton contributions have been determined in 
Ref. [12] and lead to the additional contribution c)^ (M^, Mr) 

cf{M,, Mr) = cf{M,) + cf\M,) + cf{M,) + cf{M,, M.) . (15) 

Similarly also at three-loop order the new contributions with heavy fermion loop insertions 
arise 



f (M^, Mr) = cf{M,) + cf\M,) + cf{M,) + cf\M,, Mr) 
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+ cjf (M^) + 4Um,) + cS(M^, M.). (16) 

The coefficient functions c|f''(M^, M^), cf^l{M^, M^) and Cpl{M^, M^-) will be presented 
in Section H] and have been determined with the help of the decoupling relations for the 
fine-structure constant and the fermion mass 

aK-i)(^) = e^{f,,a^-f\f^),rn,) a^-f\fi), (17) 

The decoupling functions (^{fj., a^'^f'' {ix),mh) and Cm(/^, a*^"^'^ {fj.),mh) can be obtained from 
the known QCD results ^62n69j and are given in Appendix[B]for completeness. The mass 
rrih is the mass of the decoupled heavy fermion and the superscript Uf denotes the number 
of active lepton flavors. As a check of this procedure we reproduced in Eq. f[T^ the known 
two-loop resuh of Ref. pj for c);'{M^,Mr). It can also be obtained from the results in 
Refs. [55| l60] , which incorporate the complete fermion mass dependence. In the same way 
we derived also % (M^,Mt-) and compared it with the analytic results of Ref. [60]. The 
remaining coefficient functions in Eq. ( IT6|) are again known and given in Appendix O 

In the case of the mass conversion relation from the MS to the on-shell scheme for the 
electron the internal loops in the Feynman diagrams cannot have a lower mass than the 
external fermion, since only internal electrons or heavy muons and r-leptons can appear. 
These situations are shown in Fig. HI 




Figure 4: Lepton self-energy example- diagrams at two-loop and three-loop order in QED 
which contribute to Eq. (fT2|l . The wavy lines are photons, the dashed and solid lines are 
leptons. Leptons with a dashed line have a smaller mass than leptons with a solid line. 

Similarly to Eqs. ( TT3l) and ( |T5l) the two-loop contribution reads 

42)(Me, M^, M.) = cJ^M,) + cg^Me) + (Me, M^) + (Me, M.) (19) 

and the three-loop term follows in analogy to Eqs. f[T^ and flTB]) 

cf )(Me, M^, M.) = cf{M,) + 4\m,) + (Me, M^) + cf (Me, Mr) 
+ cgl(Me, M^) + cg;(Me, Mr) + eg) (Me, M^, Mr) 
+ c%{M,) + cfl (Me, M^) + cS (Me, M.) . (20) 

The new contributions with heavy fermion loop insertions are the terms with the coefficient 
functions which have a subscript h as well as the term cj^}{Me., M^, Mr). The latter arises 
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from diagrams like the one shown in Fig. Hl^c). These contributions have again been 
obtained with the decouphng relations of Eqs. (fT7|) and ( IT8|) and are given in the next 
section. 

We have checked that the dependence of the conversion relations on the renormaliza- 
tion scale /i obeys the mass anomalous dimension [5^1701472] . 

As last ingredient the conversion factor Caa which is mediating between the fine- 
structure constant renormalized in the MS and the on-shell scheme, a{^) = CaaOt, is 
needed to order a^. It is known in literature [23l|2llllZ] to four- loop order and can 
be obtained with the help of the concept of the invariant charge [73l|71] by computing 
the vacuum polarization function at = in the MS scheme. The conversion factor 
respecting the mass hierarchy of the three charged leptons reads 



1 + 
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The symbol L^m denotes the logarithm L^m = log(/i^/M^). 

We have checked by combining the result of the vacuum polarization function in the MS 
scheme up to four-loop order at = with the above MS to on-shell conversion relations 
for the lepton masses and the fine-structure constant that Ilos{q^ = 0) is zero. We have 
also checked that the renormalization scale dependence /i of the vacuum polarization 
function cancels after the conversion to the on-shell scheme. 



4 Results 

In this section we summarize our results obtained with the methods described in the 
previous section. We begin with the terms due to heavier leptons in the conversion relation 
of the lepton masses between the MS and the on-shell scheme at three-loop order. The 
results which are known in literature are listed for completeness in Appendix O The three 
terms which contribute in Eqs. (I16p and (1201) read 
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Here we defined a„ = Li„(l) with the polylogarithm function 

X^fcli f^- The symbol is the Riemann zeta function (n = YlT=i The 
lepton mass in Eqs. (1221) - (12^ is always larger than the lepton mass M^. 

With these results one can derive the leptonic contributions to the on-shell vacuum 
polarization function from all three charged lepton generations at four-loop order in QED 
in an expansion for large external momentum. As a check of our computer routines 
which perform the algebraic manipulations we reproduce the known one-loop and two- 
loop results as well as the three-loop results U^^'{q'^,Mi), uf\q'^ , Mi, Mj) , UP{q'^,Mi), 
nf^(g2,Mi,Mj) which are given in Eqs. (5)-(10) of Ref. [12] and list them for completeness 
in Appendix |Al The new four-loop results in Eq. (fTOj) read 
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with the logarithm LqM = log(— g^/M^). The dots stand for a deeper expansion in the 
ratio of the lepton masses and the external momentum squared, M^/g^. The remaining 
symbols are the same as those defined after Eqs. fl2^ - fl2^ . 

In the following we use Eq. to perform a numerical evaluation of Aaiep by using 
the available results of the vacuum polarization function in the on-shell scheme up to 
four-loop order in order to study the relative size of the different contributions. For the 
numerical evaluation we use the following input values from CODATA [1] and PDG [75] 
for the lepton masses Me = 0.510998928(11) MeV, = 105.6583715(35) MeV, = 
1776.82(16) MeV, and for the Z-boson mass Mz = 91.1876(21) GeV. With these values 
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the coefficients of the different orders in perturbation theory of Aaiep read 

Aaiep = 13.52631(8) + (^^)^ 14.38553(6) 

+ i^-j 84.8285(7) + (^-j [810.65(1)ns - 39.8893(5)si] + 0{a^) , (35) 

where we have decomposed the new four-loop term into the contribution from non- 
singlet (ns) and singlet (SI) type of diagrams. The size of the contribution from the singlet 
diagrams is much smaller than the one from non-singlet diagrams. For the numerical 
evaluations in this section we use for the one-loop contribution the unexpanded result of 
Eq. ([6]). At two- and three- loop order we take the expanded results, where we include 
as many series coefficients in the small lepton mass expansion until their contribution 
becomes negligible compared to the numerical values given here. The same holds for the 
corresponding values in Table [H Some of the lowest expansion coefficients are shown in 
Eq. ([8]) for the two-loop order and in Eqs. fl37|) - fH0|) for the three- loop contribution. At 
four loops we use the new result up to the constant order in the mass expansion which 
was presented in this section in combination with the known terms which are summarized 
in Appendix IA.2I 

Evaluating also numerically the fine-structure constant a = 7.2973525698(24) x 10^^, 
which is taken from CODATA [1], one obtains 

Aaiep = 0.0314192... + 0.00007762... + 0.00000106... + 0.00000002... , (36) 

where each term stands for the next order in the perturbative expansion of QED, i.e. the 
one-loop, two-loop, three-loop and four-loop contribution. 

The power suppressed terms (proportional to M^/M| and higher powers) at two-loop 
order amount to about ~ 5 x 10"*^ and are smaller than the three- loop result up to the 
constant order in the small mass expansion. The latter amounts to about ~ 10^^. The 
power suppressed terms at three-loop order amount to about ~ 4 x 10~^° and are again 
smaller than the four-loop result up to the constant order, which is ~ 2 x 10~^. 

A more detailed decomposition of the individual contributions which arise from differ- 
ent gauge invariant subsets of diagrams is listed in Table [T]for the non-singlet contributions 
and in Table H] for the four- loop singlet contributions. 

At three-loop and also at four-loop order the dominant contribution to Aaiep emerges 
from diagrams with internal lepton loop insertions, especially from the insertions of elec- 
trons, which can be seen in Table [H The reason for this dominance is an logarithmic 
enhancement of these terms which arises from large logarithms in the small mass ratio 
Mg/M|. In comparison the contributions from quenched diagrams are small. This be- 
havior was already observed at three-loop order in Ref. [12] and also holds at four-loop 
order. 

Summing up all contributions one obtains Aaiep = 0.0314979 with an uncertainty of 
about ~ 2 ■ 10"'', mainly due to the error in the r-lepton mass and the Z-boson mass. 
This uncertainty is about one order of magnitude bigger than the size of the new four- 
loop contribution, which is about ~ 2 ■ 10"*^, so that the theory uncertainty in Aaiep due 
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174730.4(4) 
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Table 1: The table shows the numerical evaluation of the different gauge invariant subsets 
which contribute to Aaiep from one-loop to four-loop order. Some example-diagrams 
which correspond to the different contributions are given in Figs. [T], [2] and [31 Starting from 
three-loop order the table shows the different contributions which arise from diagrams with 
external (outer) leptons and internal (inner) leptons separately. At four-loop order arise for 
the first time diagrams with two internal lepton insertions. The non-singlet contributions 
at four-loop order are not given here, but are shown in Table [2l 
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ons in right 1 
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leptons in 
left loop 
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-0.0030194 
-0.0016177 
-0.0008029 


-0.0016177 
-0.0014803 
-0.0008029 


-0.0008029 
-0.0008029 
-0.0006656 



Table 2: Numerical evaluation of the non-singlet terms at four- loop order in QED. The 
last two diagrams of Fig. [3] are typical examples which contribute to this table. 



to the truncation of the perturbative series has been completely removed. The analytic 
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formulae which were presented in this work can be easily used to reevaluate Aaiep, once 
more precise values of the input parameters become available. 

5 Summary and conclusion 

We have determined the four-loop QED corrections to the vacuum polarization function 
in the on-shell scheme in the high energy expansion taking into account the mass hierarchy 
of all charged leptons in the Standard Model. In order to derive this on-shell result we 
combine known MS results for the vacuum polarization function in the massless limit with 
newly computed contributions at = 0, which allow to do the conversion from the MS 
to the on-shell scheme. The on-shell result is then used to perform a determination of the 
leptonic contributions to Aa which describes the running of the effective electromagnetic 
coupling to the Z-boson mass scale. The quantity Aaiep is an important ingredient for 
constraining the Standard Model in the context of electroweak precision measurements. 
We perform an numerical evaluation of Aaiep in order to study the relative size of the 
different contributions. The size of the new four-loop contribution is small and completely 
removes the theory uncertainty due to the truncation of the perturbative expansion. The 
QED corrections to Aaiep are thus from theory side well under control compared to the 
error in Aaiep arising from the uncertainty in the mass of the r-lepton and the Z-boson 
mass. 

As a by-product of the calculation we also determine the contributions due to heavy 
fermions to the conversion relation between the MS and the on-shell scheme for the lepton 
masses in QED. 
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A The on-shell vacuum polarization function in the 
high-energy limit at three-loop order and known 
four-loop results 

In this appendix we provide known results for the vacuum polarization function at three- 
loop [12j and four-loop [23l[2^ order from literature in order to have all components in a 
uniform notation. 
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A.l Three-loop results 

The results of Ref. [12] are given here for completeness and read 
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For the contributions Il^2\q'^, Mi) and Ilp\q'^, Mi) in Eqs. (1371) and (l38ll we have taken the 
sub-leading terms in the small mass expansion from Refs. [16l[l7l[19l[2T], which we have 
converted in addition from the MS to the on-shell scheme. Using q2e and exp [361438] 
we have determined also the sub-leading terms for Il['^\q'^, M^, Mj) and Il^^\q'^, Mi, Mj) 
in Eqs. ( 139|) and ( HOj) . The arising integrals have been computed in two ways: on the one 
hand they were determined by applying Laporta's algorithm [39|B0] in order to reduce all 
integrals to a set of known master integrals, and on the other hand they were calculated 
by using the MATAD and MINCER [781 - I80] routines. The dots in the equations stand for 
terms of higher order in the small mass expansion. The mass Mi in Eqs. (J39ll and (HOl) is 
the mass of the lepton in the outer loop which connects to the external photons, whereas 
the mass Mj is the mass of the internal lepton loop. The definition of the remaining 
symbols which appear in Eqs. f l37|) -f H0|) can be found in Section |H 



A. 2 Four-loop results 

The results for the vacuum polarization function of Refs. [251 [2^ are given here for com- 
pleteness and read 
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73, 4/ X 5309, 8,, 20, 146 ,^/a/2\ 
The definition of tlie different symbols can be found again in Section HI 

B Decoupling functions 

Tlie decoupling functions in Eqs. (fT7|) and (fT8|) are known in QCD to high order in 
perturbation theory [621469] and can be used for the QED calculation considered here after 
inserting the proper values for the color factors. We give the results of the perturbative 
expansion up to the order needed within this work 
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with -L^m^ = log|^^j, where ra^ is the heavy particle mass renormalized in the MS 
scheme which has been integrated out in the effective field theory. The number of active 
fermions is rif = rii + with Uh = 1 and rii is the number of light fermions, which are 
considered as massless. The definition of the remaining symbols in Eqs. (145|) and ( H6|) can 
be found in Section |3] and HI 



C The relation between the MS and on-shell scheme 
for lepton masses 

The known relations up to three-loop order in QED for the conversion of the lepton 
masses from the MS to the on-shell scheme [5H460] are given here for completeness. The 
expansion coefficient in Eq. f|T2|) at one-loop order reads 

cf\M,) = -l-h^M,. (47) 
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The coefficients in Eqs. ( !T3|l . (fT5|) and (fT9l) at two- loop order are given by 
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and at three- loop order the coefficients of Eqs. ([H]), (ITBIl and (|2U]) read 
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with Lf^M = log and > Mi. The definition of the other symbols can be found 

in Section |H The remaining three-loop coefficients which are needed for Eqs. ( IT6|) and 
(!20|) are given in Eqs. (122|) to ( 125|) of Section HI The complete mass dependence of the 
contributions from diagrams with massive lepton loop insertions at two- and three-loop 
order has been determined in Refs. [551160] . We restrict ourselves in this section to the 
contributions needed for this work. 
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